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Abstract 

In the present work, we study the cosmological model with fermionic field and with the 
non-canonical kinetic term (fermionic k-essence or /- essence). We also present some important 
reductions of the model as well as its some generalizations. We also found the exact solution 
of the model and examine the influence of such gravity-fermion interaction on the observed 
accelerated expansion of our universe. Several fermionic and bosonic-fermionic DBI models 
were constructed as some particular reductions of f-essence and g-essence, respectively and 
their Chaplygin gas counterparts are found. Finally, some nonlinear models of bouncing and 
cyclic universes are proposed. 
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1 Introduction 



The observational evidence from different sources for the present stage of accelerated expansion 
of our universe has driven the quest for theoretical explanations of such feature. Assuming the 
validity of the theory of gravity, one attempt of explanation is the existence of an unregarded, but 
dominated at present time, ingredient of the energy content of the universe, known as dark energy, 
with unusual physical properties. The other possibility is modifying the general theory of relativity 
at large scales. In cosmology, the investigation for the constituents responsible for the accelerated 
periods in the evolution of the universe is of great interest. The mysterious dark energy has been 
proposed as a cause for the late time dynamics of the current accelerated phase of the universe. 

During last years theories described by the action with non-standard kinetic terms, k-essence, 
attracted a considerable interest. Such theories were first studied in the context of k- inflation [1 , 
and then the k-essence models were suggested as dynamical dark energy for solving the cosmic 
coincidence problem [2]- [5]. The action of the k-essence scalar field <f> minimally coupled to the 
gravitational field we write in the form (see e.g. [I]-[l]) 

S = J d 4 X yf^[R + Kx{X, (/>)], (1.1) 

where 

X = O.S^V^V^, (1.2) 
is the canonical kinetic term, is the covariant derivative associated with metric Qy,w 

The 

important particular reductions of the scalar k-essence (1.1) are: 

i) K\ — Ai(X) (purely kinetic case), 

ii) K l = A 1 {X)B x {<t>), 

iii) K 1 = A 1 (X) + B 1 (4 > ). 

Note that for the FRW metric, the equations of the k-essence (1.1) have the form 

3H 2 -0.5p k = 0, (1.3) 

2H + 3H 2 + 0.5 Pk = 0, (1.4) 

K 1X 4> + (K 1X + 3HK lx )j> - K 14> = 0, (1.5) 

Pk + 'AH{ Pk + Pk ) = 0. (1.6) 

Here the kinetic term, the energy density and the pressure are given by 

X = O.50 2 , p k =2K lx X-K 1 , Vk = K x . (1.7) 

The system (1.3)-(1.6) was studied from different points of view. As one of example, we can 
consider some integrable reductions of the system (1.3)-(1.6) (see e.g. [3T]-[31]). As an example, 
let us consider the case when K\ has the form 

K x = F(X) - V{4>) (1.8) 

and we assume that F obeys one of Painleve equations e.g. the Pin - equation 

Fxx = F- X F X - X-\F X - aF 2 - j3) + 7F 3 + SF-\ (1.9) 

Then the system (1.3)-(1.6) takes the form 

3H 2 -0.5p k = 0, (1.10) 

2H + 3H 2 + 0.5p k = 0, (1.11) 

F x 4+(F x +3HFx)j> + V4, = 0, (1.12) 

F xx =F~ 1 F% -X~ l {F x - aF 2 - p)+jF 3 + 5F- 1 = 0, (1.13) 

Pk + 3H{p k + Pk ) = 0. (1.14) 
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Another class integrable reductions of the system (1.3)-(1.6) we can construct e.g. if we demand 
that <fi satisfies one of Painleve equations e.g. the Pj - equation. Then the system (1.3)-(1.6) 
becomes 

3H 2 -0.5p k = 0, (1.15) 

2H + 3H' 2 + 0.5p k = 0, (1.16) 

K lx 4>+(k lx +3HK lx )<P-K 1(t> = 0, (1.17) 

<j>HH~&4> 2 -H = 0, (1.18) 

Pk+3H( Pk + Pk ) = 0. (1.19) 

In the recent years several approaches were made to explain the accelerated expansion by choosing 
fermionic fields as the gravitational sources of energy (see e.g. refs. [6]-|25]). In particular, 
it was shown that the fermionic field plays very important role in: i) isotropization of initially 
anisotropic spacetime; ii) formation of singularity free cosmological solutions; iii) explaining late- 
time acceleration. In the present work, we study the cosmological model with fermionic field, the 
M33 - model, which has the non-canonical kinetic term ( f-essence). We examine the influence of 
such gravity-fermionic interaction on the accelerated expansion of the Universe. The formulation 
of the gravity-fermionic theory has been discussed in detail elsewhere [26]- [29]., so we will only 
present the result here. 

2 Einstein-Dirac equations 

In order to have this work self-consistent, in this section we present briefly the techniques that 
are used to include fermionic sources in the Einstein theory of gravitation and for a more detailed 
analysis the reader is referred to [5B]-[2S]. The general covariance principle imposes that the 
Dirac-Pauli matrices 7° must be replaced by their generalized counterparts T M = e^7 a , whereas 
the generalized Dirac-Pauli matrices satisfy the extended Clifford algebra, i.e., {T ll ,r 1 '} = 2g^ v . 
The Einstein-Dirac action reads as 

r J4 



S = J d 4 xy^j[R + eY - V], (2.1) 

where e = ±1 (e = 1 is the usual case and e = — 1 is the phantom case) and 

Y = 0.5i[^D^-{D^)T^}, V = V$,i;). (2.2) 

Here ip and ip = ip + ~f° denote the fermionic field and its adjoint, respectively and R is the Ricci 
scalar. The covariant derivatives are given by 

D^ = a^-%ij, D$ = d$ + $%, (2.3) 

where the spin connection is given by 

fi M = -0.25g^[r^ - ^(cU^bV, (2.4) 

with T v aX denoting the Christoffel symbols. The closed system of the equations for the model (2.1) 
looks like 



PS 



PS = V i PS 



T 


= 0, 


(2.5) 


dV 

dip 


= 0, 


(2.6) 


dV 
~chp 


= 0, 


(2.7) 


ps) 


= 0, 


(2.8) 


by 






= eY 


- V. 


(2.9) 
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3 F-essence 

Let us now we consider the M 33 - model or f-essence, which has the action 

S = j d A x^[R + K 2 {Y^M, (3.1) 
where K 2 is some function of its arguments and the canonical kinetic term has the form 

Y = 0.5i$r"D M V - (£VW r "# ( 3 - 2 ) 

We work with a space-time metric of the form 

ds 2 = -dt 2 + a 2 (dx 2 + dy 2 + dz 2 ), (3.3) 
that is the FRW metric. For this metric, the vierbein is chosen to be 

{e£) = diag(l,l/a,l/a,l/a), (eJJ) = diag(l, a, a, a). (3.4) 
The Dirac matrices of curved spacetimc T M are 

r° = 7 °, I^a-y, r 2 =a-y, r 3 = a-y, (3.5) 

r = 7°, r 1 = a 7 1 , r 2 = a 72 , r 3 -a 73 . (3.6) 

Hence we get 

O = 0, Q 1 =0.5a-y 1 -f°, fl 2 = 0.5a 7 2 7 °, 3 = 0.5a 7 3 7 °. (3.7) 
We now ready to write the equations of the M 33 - model (3.1). We have 

SH 2 + 0.5[K 2 + 0.5(K 2 ^ + K 2 ^)} = 0, (3.8) 

2H + 3H 2 + 0.5K 2 = 0, (3.9) 

K 2 Yip + 0.5(3HK 2 Y + K2Y)i>-h K 2 j = 0, (3.10) 

K 2Y tP + 0.5(3HK 2Y + K 2Y )i> + iK 2i)1 = 0, (3.11) 

Pf + mpf+Pf) = 0, (3.12) 

where Y = 0.5z(i/> 7 ( V — i/) 7 ^) and 

p f =K 2Y Y-K 2 = -[K 2 + 0.5{K 2 ^ + K 2 ^)}, p f = K 2 (3.13) 

are the energy density and pressure of the fermionic field. If K 2 — Y — V, then from the system 
(3.8)-(3.12) we get the corresponding equations of the Einstcin-Dirac model. 

3.1 Submodels 

3.1.1 The M 33A - model 

The Langrangian of the M 33 a - model has the form 

K 2 = A 2 (Y). (3.14) 
So the M 33j 4 - model is the purely kinetic fermionic k-essence. 

3.1.2 The M 33B - model 

The Langrangian of the M 33B - model reads as 

K 2 = A 2 {Y)B 2 {i>,i,). (3.15) 
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3.1.3 The M 33C - model 

The Langrangian of the M 33 a - model has the form 

K 2 = A 2 (Y) + B 2 (iP,ijj). 

3.2 Solution 



(3.16) 



In this subsection we want to construct a solution of the M33-model. Let K 2 = K 2 (Y,u), where 
u = tp^p. Then the system (3.8)-(3.12) becomes 



3H 2 + 0.5[K 2 +K' 2 u] 


= o, 


(3.17) 


2H + 3H 2 + 0.5K 2 


= o, 


(3.18) 


K 2Y ip + 0.5(3HK 2Y + K 2Y )^ - h°K 2 iP 


= o, 


(3.19) 


K 2Y $ + 0.5(3HK 2Y + K 2Y )4> + iK 2 i)f 


= o, 


(3.20) 


p f + 3H( Pf + Pf ) 


= o, 


(3.21) 



where K 2 = dK 2 /du and 



Pf = - [K 2 + K 2 u] . p f = 



(3.22) 



We now consider the submodel (3.16) that is the M 3 3c-model, where we assume that A 2 = 
aY n , B 2 = (3u m that is the case K 2 = aY n + j3u m . Let a = a t x . Then we have the following 
solution 



Y 



6mA 2 -4(m + l)A 



am 



tr 



»1 



where 



A = 



2n — 2m + 2mn 
3mn 



3/ c 

V an 



6mA 2 — 4(m - 


-1)A" 




4A" 


am 







Finally we present the following formulas 



(3.23) 



(3.24) 



(3.25) 



where c = \c\\ + \c\ 2 — | c|§ — \c\\, Cj = consts. 
3.3 f-DBI 

Let us rewrite the action (3.1) as 

S = J d i x^[R + 2K 2 {Y^,^)]. 

For the FRW metric the corresponding system of equations takes the form 

3H 2 - Pf = 0, 
2H + 3H 2 +p f = 0, 
K 2Y i) + §.h(3HK 2Y + k 2Y )^-i^K 2 ^ = 0, 

K 2Y iP + 0.5(3HK 2Y + k 2Y )i> + iK 2 ^ = 0, 
Pf + 3H(pf+ Pf ) = 0, 

where Y = 0.5i(i/>7 C V — Tp^ip) and 

p f = YK 2Y - K 2 , pf = K 2 . 
Now we want present some examples fermionic DBI models or short f-DBI. 



(3.26) 



(3.27) 
(3.28) 
(3.29) 

(3.30) 
(3.31) 



(3.32) 
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3.3.1 Tachyonic models 

i) First let us consider the following tachyonic models 



K 2 = -tryi - aY, 

where 

U=U(f,T), Y = 0.5i(fj°f -T^T) 
and T is the fermionic tachyon field. Then we get the following EoS 



P = -P T V P 2 - U 2 , 

where 

0.5U 



VI -aY 

ii) Let the f-essence Lagrangian has the form 



+ O.bUVl - aY. 



K 2 = -U\J\ - aY 2 , 

where 

U=U(f,T), Y = 0.5i(fj°f-f 1 °T) 
and T is the fermionic tachyon field. Then we get the following EoS 

U 2 

P = , 

P 

where 

U 



Vl - aY 2 

So this model corresponds to the Chaplygin gas [3"U] . 

3.3.2 Dark energy models 

a) Now we want consider the following DBI model 



K 3 = eU(^l-a^-l) + V, 
where 

U = U($,ip), Y = 0.5i(^7°V) - ^7V), V = V($,ip) 
and t/j is the classical commuting fermionic field. Then we get the following EoS 



p = -2p ± VV - 2e 2 t/ 2 + F - eJJ, 
where p — p — ell + V and the energy density is given by 



y n ^fU 
p = -0.5eU\l I- a + eU —V. 

1 



b) Our second example reads as 



K 2 = eU(<Jl-a^--l) + V, 
where 

U = U(ip,ip), Y = 0.5i(V>7°^ — ^j°ip), V = V(i>,ip) 
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and ip is the classical commuting fermionic field. Then we get the following EoS 



P=-J^y-*> + V, (3.47) 



where the energy density is given by 



p = ; eU +eU -V. (3.48) 

1 Y 2 

1-a— 



From (3.47) follows that as V — eU this model becomes 

e 2 U 2 

p = — (3.49) 

P 

that corresponds to the Chaplygin gas [30] . 
3.4 Integrable f-essence models 

One of interesting class of f-essence models is integrable models (see e.g. [3T]-[33]). Let the 
Lagrangian of f-essence has the form 

K 2 =F(Y)-V(ij>,ip). (3.50) 

As an example, here we consider the following particular f-essence model when F(Y) obeys the 
equation 

F YY = 2F 3 + YF + a, (3.51) 

where a is some constant. It is nothing but the P// - equation so that is integrable. Now the 
equations off-essence (3.27)-(3.31) take the form 

(3.52) 
(3.53) 
(3.54) 

(3.55) 
(3.56) 
(3.57) 

To construct solutions of this system we start from the equation (3.56). The equation (3.56) has 
the following particular solutions (see e.g. |31)-|34]) 

F = F(Y;l.5) = ip-(2i; 2 +Y)-\ (3.58) 

F ee F(y;l)=-i (3.59) 
1 w 2 

F = F ^ Y ^=Y-WT-4' (3 - 60) 

F F(Y- 3)- 6y2 ( r3 + 1 °) refill 

= F(y ' 3) -y^TI-r6 + 20y 3_ 80 ' ( 3 - 61 ) 

F F(Y-4) 1 I 6y2 ( y3 + 1 °) 9r 5 (y 3 +40) 

[ ' > y r 6 + 2oy 3 -80 r 9 + 6oy 6 + 11200' K ' 

F ee F(Y; 0.5e) = (3.63) 

and so on. Here 

ip = (ln0)y, <j)(Y) = Ci^(-2- 1 / 3 r) + C 2 ^(-2- 1 / 3 F), (3.64) 
where C% — consts and Ai(x), Bi(x) are Airy functions. 





3H 2 - Pf 


= o, 




2H + 3H 2 + Pf 


= o, 


F Y ip + 0.5(3HF Y -t 


-F Y )Tp + ij% 


- o, 


F Y 4> + 0.5(3HF Y -t 


- F Y )i>-iV^° 


= o, 


Fyy — 


2F 3 -YF -a 


= o, 


Pf 


+ 3H(p f + Pf ) 


= 0. 
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4 G-essence 

4.1 Equations 

We now would like to present the M34 - model (the generalized k-essence or g-essence for short) 
which has the following action 

S = j d A x^[R + 2K{X,Y^^M- (4-1) 

For the FRW metric (3.3), the equations of g-essence (4.1) have the form 

3H 2 -p = 0, (4.2) 

2H + 3H 2 +p = 0, (4.3) 

Kx$+{Kx + 3HK x )4>-K4, = 0, (4.4) 

KYip + O.biiHKY + Ky^-H ^ = 0, (4.5) 

Kyip + 0.5('SHKy + Ky)iP + iK^° = 0, (4.6) 

p + 3H(p + p) = 0. (4.7) 

Here the kinetic terms, the energy density and the pressure are given by 

X = 0.5<j> 2 , Y = 0.5i(^y°V> ~ V¥V) (4-8) 

and 

p = 2XK X + YK Y -K, p = K. (4.9) 

Note that the model (4.1) contents some important particular submodels. For example: 

i) the scalar k-essence (1.1) as K = Ki(X,<j>); 

ii) the M 33 - model (3.1) as K = K 2 (Y,ip,^); 

hi) the M 3iA - model as K = K X {X, <f>)K 2 (Y, ip, ip);_ 
vi) the M 34i j - model as K = K X {X, <j>) + K 2 (Y, ip, # 

4.2 Bosonic-fermionic Chaplygin gas 

As for k-essence and f-essence cases we can find the bosonic-fermionic Chaplygin gas from g-essence. 
To do it, let us consider the following particular g-essence model 

K = Vb[l + VxX 1 ' 271 + V2Y 1 ' 71 + v 3 X v Y^- 2nv ^ n ] n , (4.10) 

where Vj = Vj(<p, ip, tp) and n, v = consts. Then we have 

p = Voil + ViX 1 / 2 * 1 + V2Y 1 ' 11 + V s X v Y^- 2n ^' n ] n , (4.11) 

p = -Voll + VxX 1 ' 271 + V 2 Y 1 ' n + v 3 X u Y {l - 2nu ^ n ] n - 1 . (4.12) 

Hence we get the generalized bosonic-fermionic Chaplygin gas 

pl-n p a 

where a = n/(l — n), A = (— Vo) 1 ^ . The case n — 0.5 corresponds to the bosonic-fermionic 
Chaplygin gas. 

4.3 Bosonic-fermionic DBI 

Consider the following bosonic-fermionic DBI model 

K = V [l + yiX 1/2n + V 2 Y 1/n + v 3 X u Y^- 2nu ^ n ] n + Va, (4.14) 
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where Vj = Vj((f>, if>, ip) and n, v — consts. If n = 0.5 then it takes the form 



K = Wl + ViX + V 2 Y 2 + VzX»Y 2 ( l -») + V 4 . (4.15) 

From (4.14) we have 

p = V [l + V!X 1/2n + V 2 Y 1/n + y 3 x u Y {1 - 2nu)/n ] n + V 4 , (4.16) 

p = -Vail + VxX 1 ' 2 ™ + V 2 Y 1 ' n + V z X»Y {1 - 2nv V n ] n - 1 -Vi. (4.17) 

Hence we get the generalized bosonic-fermionic Chaplygin gas 

where a = n/(l — n), A = (— Vo) ■ The case n = 0.5 corresponds to the bosonic-fermionic 
Chaplygin gas. The interesting particular case of the model (4.14) we get if we set 

K = V {[1 + Vi'V^X 1 / 2 " + V^VqY 1 ' 71 + v' z v' Q X v Y {1 ~ 2n ^l n \ n - 1} + V(&&^), (4.19) 

where V^- = (0, ?/>, V 7 )- Let n = 0.5. Then the bosonic-fermionic DBI (4.18) becomes 



K = Vo[Vl + PiX + foY 2 + p 3 X"Y^- v ) - 1] + V(cf>, i/>, (4.20) 
where /3 3 = V-V^ 1 . 

5 On some models of bouncing and cyclic universes 

Finally we would like to note that it is interesting to find some f-essence models which describe 
the cyclic and bouncing universes. Here we briefly present some of such models (see e.g. |3T ) -[3"I | ). 
i) Consider the following model 

at — asggt + 2k<2£ + 3aa£ — 2a^a^ — aa^ = (5-1) 

or its twin 

ei£ - a«£ + 2Ka t + 3aa t - 2a t a tt - aa ttt = (5.2) 

which is nothing but the Camassa-Holm equation. Here a = a(t, £) is the scale factor, k = const and 
£ is some artificial coordinate e.g. £ = A (cosmological constant), £ = G (Newton's gravitational 
constant), £, — k (spatial curvature) or £ = c (speed of light) see e.g. 

a 2 8ttG kc 2 Ac 2 a AkG , 9 , Ac 2 

- =^P-- + X' - = ^-(P + 3c- 2 p) + — • (5.3) 

a 3 or 3 a 3 3 

In the case of no linear dispersion, k = and £ = A, the equations (5.1) and (5.2) possess singular 
solutions in the form of peaked travelling waves called peakons 

a = 0.5f]e^-^ (5.4) 

and 



a = O^e-I*-^ 1 , (5.5) 

respectively. They are 1-peakon solutions. The multi-peakon solution versions of these solutions 
have the form 

n 

a = ^2r ] {t)e- 2 ^- q ^ (5.6) 



3=1 

and 



2>(0e- a l*-*<0l, (5.7) 
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respectively. Our idea is that the n-peakon solutions can describe the cyclic and bounce universes. 
Hope that we will return to this question in the future from the more detail position. 

ii) Our second example is the following model [here we set e.g. £ = A and a a = da/dA, a t = 
da/dt etc]: 

a t = 0.75(a 2 ) A (5.8) 

or its twin 

a A = 0.75(a 2 ) t (5.9) 

which are known to develop shocks. Eqs. (5.8) and (5.9) are nothing but the dispersionless 
Korteweg-de Vries equations (dKdVE) or Riemann equations. It is well-known that the equations 
(5.8) and (5.9) have the following solutions 

<z(A,t) = /i(A-0.75<rf) (5.10) 

and 

a(t,A) = h(t-0.75aA), (5.11) 

respectively. Here h is an arbitrary function. Also we see e.g. from the solution (5.10) that the 
velocity of a point of the wave, with constant amplitude a, is proportional to its amplitude leading 
to the "breaking" of the wave. Also we note that the wave also develops discontinuities in its 
evolution. Let's now we give the following particular solutions of the equations (5.8) and (5.9): 

a=(/3 1 +/3 2 A)(/3 3 -1.5/3 2 i)- 1 (5.12) 

and 

a = (Px + p 2 t)(p 3 ~ 1.5/3 2 A)-\ (5.13) 

respectively, where /3, = const s. 



6 Conclusion 

We briefly summarize the present work. We first derived the equations of the M33 - model for the 
FRW space-time. Then we found their exact solution for the /-essence AT 2 = aY n + (3u m . Finally, 
let us we present the expressions for the equation of state and deceleration parameters. For the 
our particular solution (3.19) they take the form 

mn 3mn 

Wf =-l-\ , q= -I + 777 7 ■ (6.1) 

n — ui -T mn 2(n — m + mn) 

These formulas tell us that the M33 - model can describes the observed accelerated expansion of 
our universe. Also we have been constructed the fermionic DBI models from f-essence. Finally, 
some nonlinear models of bouncing and cyclic universes are proposed. 
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